EXACTLY SOLVABLE INTERACTING SPIN-ICE VERTEX MODEL 
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A special family of solvable five-vertex model is introduced on a square lattice. In addition to the 
usual nearest neighbor interactions, the vertices defining the model also interact along one of the 
diagonals of the lattice. Such family of models includes in a special limit the standard six-vertex 
model. The exact solution of these models gives the first application of the matrix product ansatz 
introduced recently and applied successfully in the solution of quantum chains. The phase diagram 
and the free energy of the models are calculated in the thermodynamic limit. The models exhibit 
massless phases and our analytical and numerical analysis indicate that such phases are governed 
by a conformal field theory with central charge c = 1 and continuosly varying critical exponents. 
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I. INTRODUCTION 

The six- vertex model was introduced by Pauling in or- 
der to explain the residual entropy of ice at zero temper- 
ature. Due to its exact integrability and its connection 
with the XXZ quantum chain, this vertex model is con- 
sidered as a paradigm of exact integrability 0, 0, . In 
this model the interactions are between nearest neighbor 
vertices and are ruled by the geometrical connectivity of 
the allowed vertex configurations on the lattice. The in- 
teraction energy among two vertices is zero for allowed 
configurations and infinite otherwise. In this paper we 
introduce a special family of five-vertex models. Besides 
the usual nearest-neighbor interactions, imposed by their 
connectivity, there exist additional interactions among 
vertices at larger distances. Moreover we are going to 
show that this family of models is exactly integrable and 
contains as a special case the standard six-vertex model. 

Although the introduction of a new exactly integrable 
model is interesting by itself we would like to mention 
that recently the interest on two-dimensional vertex mod- 
els was rised by their possible physical realizations 0. 
Artificially produced materials have exhibited a frus- 
trated square-lattice array of nanosize ferromagnetic do- 
mains. These materials, for large lattice spacing, are in 
an ordered "spin- ice" square lattice where their bonds 
satisfy the six vertex connectivity rules. By decreasing 
the lattice spacing, the increased frustration produces ex- 
tra interactions among the vertices as well as other vertex 
configurations |^|. 

The solution of the models presented on this paper 
are obtained through the exact diagonalization of the 
diagonal-to-diagonal transfer matrix. The exact solu- 
tion of transfer matrices associated to vertex models or 
quantum Hamiltonians are usually obtained through the 
Bethe ansatz |5j on its several formulations. That ansatz 
asserts that the amplitudes of the eigenfunctions of these 
operators are given by a sum of appropriate plane waves. 



Instead of making use of the Bethe ansatz, the solu- 
tion we derive will be obtained through the matrix prod- 
uct ansatz recently introduced in [(|. According to this 
ansatz, the amplitudes of the eigenfunctions are given in 
terms of a matrix product whose matrices obey special 
algebraic relations. The present paper presents the first 
application of the matrix product ansatz for the exact 
solution of a transfer matrix. 

The layout of this papers is as follows. In Sections II 
and III we introduce the solvable family of interacting 
vertex models and their associated diagonal-to-diagonal 
transfer matrix. In Sec. IV we present the matrix prod- 
uct ansatz and diagonalize the transfer matrix obtained 
in the previous section. Section V presents the general 
picture of the roots of the spectral parameter equations 
associated to the exact solution of the models. The free 
energy in the thermodynamic limit is obtained in Sec. VI 
whilst in Sec. VII the operator content of the model on 
its massless phase is obtained by exploring the conformal 
invariance of the infinite system. Finally in Sec. VIII we 
conclude our paper with a general discussion. 



II. THE INTERACTING FIVE- VERTEX 
MODEL 

The family of vertex models we introduce and solve are 
defined on a square lattice with M rows and L columns 
(see Fig. QJi). At each horizontal (vertical) link we at- 
tach an arrow pointing to the left or right (up or down) 
direction. Such arrow configurations can be equivalently 
described by the vertex configurations of the lattice. A 
vertex configuration at a given site (center) is formed 
by the four arrows attached to its links. Similar to the 
six-vertex model we impose that the allowed arrow con- 
figurations only contain vertices satisfying the ice rules, 
namely, the fugacity of a given vertex is infinite unless 
two of its four arrows point inward and the two oth- 
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ers point outward of its center. In Fig. [5] we show the 
six vertex configurations respecting the ice rules, with 
their respective fugacities ci, C2, 61, b%, ao, ai . The parti- 
tion function is given by the sum of all possible vertex 
configurations with the Boltzmann weights given by the 
product of the fugacities of the vertices. It is impor- 
tant to notice that the ice rules imply the existence of 
interactions among the nearest-neighbor vertices. These 
interactions have a zero or infinite value. For example 
(see Fig. |2J, the vertex with fugacity ao has an infinite 
interaction energy if the vertex on its right side is one 
of the vertices (ai, 61,02) or the vertex on its left side is 
one of the vertices (ai, b\,c\). The model with no extra 
interactions, besides those given by the connectivity of 
the vertices on the lattice, is the well known six-vertex 
model. This model is exactly integrable for arbitrary val- 
ues of the fugacities, and is considered a prototype of an 
exact integrable model 0, 0] • 
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FIG. 1: (a) Square lattice with M horizontal (L vertical) 
lines. The interacting five- vertex models have extra interac- 
tions along the dashed diagonals, (b) The deformed lattice. 
The dashed diagonals are now in the horizontal direction. 



We consider a special family of interacting five-vertex 
models. Besides the previously mentioned nearest- 
neighbor interactions (infinite or zero) imposed by the 
lattice connectivity, the models also contain interactions 
among pairs of vertices at larger distances. The allowed 
vertex configurations, with their respective configura- 
tions are the first five configurations shown in Fig. [21 
Contrary to the six-vertex model the vertex configura- 
tions with fugacity a\ are forbidden (zero fugacity). 

Such interacting five-vertex models are labelled by a 
fixed positive integer t that may take the values t = 
1,2,3,... . This parameter specifies the additional in- 
teractions among the vertices. These interactions occur 
along the diagonals of the square lattice that go from the 
top left to the botton right direction (see the dashed diag- 
onals in Fig.^,). A pair of vertices at distance D = l\/~2 
(I = 1, 2, . . .), in units of lattice spacing, along this diag- 
onal interacts as follows 

a) the interaction energy is zero if I > t, 

b) if one of the vertices is ao the interaction energy is 
zero for all values of I, 

c) if neither of the vertices is ao, the interaction energy 
is infinite if I < t, except in the special case where I = t 
and C2 is on the left of c\. In this case the interaction 
energy e/ is finite and produces a Boltzmann weight cj 
given by 



cic 2 



(i) 



In Fig. and Fig. [SJj we show, for the model with 
t = 2, some examples of forbidden (infinite energy) and 
allowed configurations, respectively. We also show, in 
those figures, the total contribution of the pair of vertices 
to the Boltzmann factor. In general, the contribution of 
a given pair of vertices is zero if the pair is not allowed 
(infinite interaction energy) or is given by the product of 
their fugacities (zero interaction energy). The exception 
to this rule happens when we have the vertices c\ and C2 
at the distance D = y2i along the diagonal, with C2 on 
the left of ci (see Fig. |3Jd). In this case, from Q, the 
contribution is given by C1C2C/ = a\. 

The partition function of the vertex models are usu- 
ally calculated by the diagonalization of the row-to-row 
transfer matrix, connecting the arrow configurations of 
two consecutive horizontal lines. In the case under con- 
sideration, due to the existence of the extra interactions 
along the diagonals, the calculation is much simpler if we 
consider the diagonal-to-diagonal transfer matrix, con- 
necting the arrow configurations of two consecutive diag- 
onals of the lattice. 



III. THE DIAGONAL-TO-DIAGONAL 
TRANSFER MATRIX 



FIG. 2: The six vertex configurations allowed by the ice rules, 
with their respective fugacities. 



Following Bariev [j| in order to construct the diagonal- 
to-diagonal transfer matrix for the interacting five- vertex 
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FIG. 3: Examples of forbidden (a) and allowed (b) config- 
urations of pairs of vertices for the interacting vertex model 
with t = 2. The contribution of the vertices to the Boltzmann 
weights are also shown. 



models it is convenient to distort the square lattice shown 
in Fig.QJi as in Fig.^J). In this case the vertices which are 
at closest distances along the dashed diagonals of Fig.^, 
are now at the closest distance along the horizontal di- 
rection. We are going to solve the model with toroidal 
boundary conditions on the distorted lattice of Fig. ^p. 
In the undeformed lattice of Fig. this boundary con- 
dition translates into a helical one where we identify the 
sites 



(ij) = {i + M,j) = (i-L,j + L), 



(2) 



where denote the sites located at the columns j 

(j = 1, 2, . . . , L) and rows i (i = 1, 2, . . . , M). 

In the distorted lattice of Fig. the allowed five ver- 
tex configurations of Fig. [5] are shown in Fig. 21 We also 
represent on this last figure a convenient vertex repre- 
sentation where we only draw the arrows pointing to the 
down direction. Hereafter we are going to use such rep- 
resentation. We see from Fig. 0] that the allowed vertex 
configurations have or 2 arrows. We describe a given 
configuration \ip > of n arrows on a horizontal line of the 
lattice of Fig. ^3 by given their coordinates (x\, . . . , x n ) 
and a set of labels (ax, . . . , a n ), where on — 1 or on — 2 
depending if we have at the coordinates Xi a vertical or 
inclined arrow respectively, i. e., 



\(p >= \xi,ai;x 2 ,a 2 ; 



(3) 



The diagonal-to-diagonal transfer matrix T^-d, rep- 
resented pictorically in Fig. [5] connects the arrow con- 
figurations of two horizontal lines in Fig. \\]p. Its matrix 
elements give the contribution of these configurations to 
the Boltzmann weight. For the allowed configurations 
they are given by 




FIG. 4: The five allowed vertex configurations in the dis- 
torted diagonal lattice with their respective fugacities. In the 
second line we give a representation where it is only drawn 
the arrows pointing down. 



< <p\T D - D \<f/ >= cT^b^a^c^ 



(4) 



where rii are the number of vertices of type i (i = 1, . . . , 5) 
and n p is the number of pairs of vertices c 2 — Ci at distance 
of t units along the horizontal lines (see Fig. [3J). The 
partition function is given by the trace 
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FIG. 5: Pictoric representation of the diagonal-to-diagonal 
transfer matrix. The coordinates (i, a) specify the possible 
arrows positions. 
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FIG. 6: Example of an allowed (a) and not allowed (b) con- 
figuration of the interacting five- vertex model with t = 1. The 
arrows have an effective hard-core size s — 2t + 1 = 3. 
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Z = Tr(Tjf_ D ). 



(5) 



It is interesting to observe that in the distorted lattice 
the allowed arrow configurations on the horizontal lines, 
besides having a fixed number of arrows, could also be 
interpreted as if the arrows would have an effective size 
s = 2t + 1 (s — 3, 5, . . .), in lattice spacing units of the 
distorted lattice. An arrow on a given link has hard-core 
interactions that exclude the occupation of other arrows 
at the link itself as well as at the 2i-nearest links on 
its right. In Fig. [5|we represent pictorically an allowed 
(a) and a not allowed (b) configuration of arrows for the 
interacting five- vertex model with t = 1. This example 
corresponds to arrows of effective size s = 3. 

The above interpretation, where we associate an ef- 
fective size to the arrows, allows us a simple extension 
of our model to the case where t = 0. In this case the 
arrows have a unity size and the extra hard-core interac- 
tions among the vertices with fugacities C2 and c\ occur 
when the arrows are at the same site, thus giving rise to 
an extra vertex with fugacity c\Ci^^ = ax, and hence 
the model reduces to the well known six-vertex model, 
with fugacities given in Fig. [21 



IV. THE MATRIX PRODUCT ANSATZ AND 
THE DIAGONALIZATION OF THE TRANSFER 
MATRIX T D - D 

In this section we solve the eigenvalue equation for the 
diagonal-to-diagonal transfer matrix Tu-d introduced in 
the previous section. As a consequence of the arrow con- 
servation and the translation invariance of the arrow con- 
figurations on the horizontal lines of Fig. the matrix 
Td-d can be splitted into block disjoint sectors labelled 
by the number n (n = 0, 1, ... , 2L) of arrows and mo- 
mentum p = ^jjj {j = 0, 1, . . . , L — 1). We want to solve, 
in each of these sectors, the eigenvalue equation 



An.pl^n.p >— Td- 



D\^n,p >, 



(6) 



where A„ iP and \^ n ,p > are the eigenvalues and eigen- 
vectors of Tu-D-, respectively. These eigenvectors can be 
written in general as 



> 



y^y^d aS x ^ ■ ■ -> x n)\xi,ai; ...;x n ,a n >, 



(7) 



where <\P a a ( x ii ■ ■ ■ i x n) is the amplitude correspond- 
ing to the arrow configuration where n arrows of type 
(ax, ■ ■ ■ , ot n ) are located at (xx, . . . , x n ), respectively. 
The symbol (*) in Q means that the sums of {x} and 



{a} are restricted to the sets obeying the hard-core ex- 
clusions, for a given interacting parameter t: 



Xi+i > Xi + t + 1 - S (8) 

for i = 1, . . . , n — 1, and 

t + 1 — 5 ait x5 an ,2 < x n -xx < L-t - l + 5 any iS au 2. (9) 

Since \^ n ,p > is also an eigenvector with momentum p, 
then the amplitudes satisfy 



(xi , . . . , X n ) 



Xxx + 1, ...,x n + 1) 
for x n < L — 1, whilst for x n = L 



(xx 



- ip 



i) 



(10) 



(ii) 



The exact solution of © is obtained by an appropriate 
ansatz for the unknown amplitudes <$P a a ^ (xx, • ■ ■ , x n ). 
As shown in the last section our model reduces to the 
standard six-vertex model for the case where t = (s = 
1). An appropriate coordinate Bethe ansatz, that solves 
the eigenvalue equation ©, is known in this case 0. 
Thus the amplitudes {(f%, 1 ..., a „(^i, ■ ■ ■ , x n )} are given, as 
usual, by a combination of plane waves whose wavenum- 
bers are fixed by the eigenvalue equation JBJ). 

In this paper we solve for general values of t 
(t = 0, 1, 2, . . .) or s = 2t + 1 (s = 1, 3, 5, . . .), by using a 
distinct ansatz. The matrix product ansatz we are going 
to use was introduced in [(| for exact integrable quantum 
chains. We present in this paper the first application of 
this matrix product ansatz for transfer matrices. Accord- 
ing to such ansatz the amplitudes <$? a ... a (^lj ■ ■ ■ >%n) 
are given in terms of a matrix product satisfying an un- 
known associative algebra. The model is exactly inte- 
grable if the eigenvalue equations fix consistently the al- 
gebraic relations among the matrices. 

In order to formulate the matrix product ansatz we 
make a one-to-one correspondence between configura- 
tions of arrows and products of matrices. The matrix 
product associated to a given arrow configuration is ob- 
tained by associating a matrix E to the sites with no 
arrow, a matrix to the sites with a single arrow of 
type a (a = 1, 2), and finally the matrix A^E^A^ to 
the sites with two arrows [l0| . The unknown amplitudes 
in J7J are obtained by associating them to the matrix 
product ansatz 



4%i,...,a n (*l> ■ ■ ■ > x n) & E^A^E^-^A^ 
_ -gix n -x n -i-l j^(a n ) jjjL-Xn _ (12) 

Actually E and A are abstract operators with an asso- 
ciative product. A well defined eigenfunction is obtained, 
apart from a normalization factor, if all the amplitudes 
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are uniquely related. Equivalently, in the subset of words 
(products of matrices) of the algebra containing n matri- 
ces A and (L — n) matrices E, there exists only a single 
independent word. The relation between any two words 
gives the ratio between the corresponding amplitudes in 

G2J- 

We could also formulate the ansatz l|12l) by associating 
a complex number to the single independent word. We 
can choose any operation for the matrix products that 
gives a non-zero scalar. In the original formulation of the 
matrix product ansatz for exactly integrable quantum 
chains with periodic boundary conditions [||, the trace 
operation was chosen to produce the following scalar 



, iCtn (xi, x n ) = Tr^ 1 - 1 E 3 "-" 1 - 1 

(13) 



-l A (a n ) E L- x „ n 



The matrix f2 p was chosen to have a given algebraic re- 
lation with the matrices E and A. Recently, Golinelli 
and Mallick have show that in the particular case 
of the asymmetric exclusion problem in a periodic chain 
it is possible to formulate the ansatz only by imposing 
relations between the matrices E and £l p . The relations 
between A and f2 p is then totally arbitrary. Actually, as 
we are going to show in this paper, we do not need to im- 
pose any algebraic relation between the matrices E and 
A with f2 p . The matrix fi p can be any arbitrary matrix 
that produces a non vanishing trace in (|13l) . This obser- 
vation is not particular to the present model. It is valid 
for any of the exactly integrable quantum chains solved 
in the original formulation of the matrix product ansatz 
jfj . Instead of restricting our matrix product ansatz with 
the trace operation, as in i|13fl . we will consider the more 
general formulation given by (|T2|> . 

Since the eigenfunctions produced by the ansatz have 
a well defined momentum p — ^j-j (J = 0, . . . , L— 1), the 
relations (|10(l and imply the following constraints 
for the matrix products appearing in the ansatz (|12|) 

E x l — 1 x\ — 1 i i _ ^(a„)gL- x n _ 

e ~ip-gX! ^(ai) £jX2—xi—l . , , A (a n ) jgL—x n —l ^4) 

for x n < L — 1 , and for x n = L 

JTJxi-l j^(ai) £X 2 -xi-l , _ _ A (ot„) _ 

e -ip A (a„) E xx-l A ( ai ) . _ A ( an - 1 ) E L-x n - 1 -l_ ^ 

In order to proceed, in the usual way, we are going 
to consider firstly the eigensectors of Td-d with small 
values of n. 

n = 0. We have the trivial solution presenting a single 
eigenvalue Ao = aft. 

n = l. We have in this case 2L possible configurations, 
corresponding to a single arrow of type a = 1 , 2 at any 
lattice position x = 1, . . . , L. The eigenvalue equation 
© gives us the 2L equations 



]jyc-l A (l) E L-x 



b 2 E x - 1 A^ E L ~ X +c 1 E x - 1 A^E L - x , (16) 
for x = 1, . . . , L, 

_^_ E x-l A (2) E L- X = 

c 2 E x A {1) E L - X - Y + hE x A {2) E L - x -\ (17) 
for x = 1, . . . , L — 1, and for x = L 
Ai 



_ E L-1 A (2) = C2A (D E L-1 



+hA^E L -\ 



(18) 



These equations are written in a convenient form by im- 
posing that both matrices A a (a = 1,2) depend on a 
single spectral parameter matrix A k : 



A« = t&AtE 1 -*, A^ = ^E 1 - 



(19) 



where (f>\ and 4>\ are unknown scalars and A k satisfies 
the following algebraic relation with the matrix E 



EA k = e lk A k E. 
Inserting (fH3|) on (fl?)|) we obtain 



xE x - 1 A k E 1 - 2t E L - x = 0. 



(20) 



(21) 



On the other hand inserting (|19fl on l|17|) - (|18fl and using 
(|20|l we obtain 



/ a il Z/— 1/1 ik j L — 111 ik\ 

(Ai0 2 - c 2 a <p Y e - 61 a 2 e ) 



xE x ~ 1 A k E 1 - 2t E L - x = 0, 



(22) 



for x = 1, . . . , L — 1 and 



xA k E L - 2t = 0. (23) 

However by inserting (|2t)|l in (|14|l we verify that k coin- 
cide with the momentum of the eigenstate: 



k = p 



2irm 



, m = 0,1,..., L- 1. (24) 



In order to produce a non-zero norm state we should 
impose E^ 1 A k E 1 - 2t E L ~ x ^ 0, for x = 1, . . . , L. Due 
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to the equality e %kL = 1, equations H21fl - I|23[l reduce to In the case where x 2 = L, equations (|29|l and l|30[l are 
the eigenvalue problem of a 2 x 2 matrix replaced by 



Ai 



c 2 e 



ik 



Cl 



(25) 



The diagonalization of l|25|l gives us the eigenvalue 



A 1 (fc) = A^ ) (fc) = 



L-l 



-(b 2 + b 1 e lk 



+l[(b 2 + he i y - 4e«(& 2 bi - c 2 ci)]5), (26) 
with I = ±1 and k given by 124|) . 

n = 2. The eigenvalue equation 10 gives two types 
of relations for the amplitudes l|12fl . The configuration 
where we have no collisions, i. e, \x\, a\) x 2 , a 2 
\xi, 2;Xx + t + 1, 1 > gives generalizations of the equa- 
tions (f 1 f>|> - (f 1 H|) for two particles 



E Xl ~ 1 e X2 ~ Xi ~ 1 E L ~ X2 = 



-L-2 



b\E Xl ~ x A^ e X2 ~ Xi ~ 1 A^ E L ~ X2 

+c l b 2 E Xl - 1 A {2) E X2 - Xl - 1 A {1) E L - X2 

+b 2 c 1 E Xl - 1 A ( - 1) E X2 ~ Xl - 1 A (2) E L - X2 

+ c \e Xi - 1 a (2) e X2 - Xi - 1 a {2) e l - X2 , 



A- 



j\i 2 ) JgX2-Xl—l ^(1) £L-X2 _ 



c 2 b 2 E xl A w E X2 - xl - 2 A {1) E L ~ X2 
+hb 2 E Xl A (2) E X2 - X1 - 2 A (1) E L -' X2 
+c 2Cl [E Xl A {1) E X2 - X1 - 2 A (2) E L - X2 
+61C1 E x 1 A (2) E X2 - Xl - 2 A (2) E L -' X2 , 



A. 



L-2 



c 2 b 2 E Xl - l A {1) E X2 - Xl A {1) E L - X2 - 1 
+c 1 c 2 [E Xl - 1 A {2) E X2 - Xl A {1) E L - X2 - 1 
+b 2 b 1 E X1 - 1 A (1) E X2 ~ Xl A (2) E L - X2 ~ X 
+c 1 b 1 E xl - 1 A {2) E X2 - xl A ( - 2 '>E L - X2 - 1 , 



A, 



-L-2 



^xx-lj^X) £X 2 -xi-l A (2) £ 



L—x 2 



c 2 E xi A (l) E x 2 -x 1 -l A (l) E L-x 2 -l 

+b 1 c 2 E Xl A (2) E X2 - Xl - 1 A (1) E L - X2 ~ 1 
+c 2 b 1 E Xl A {1) E X2 - Xl - 1 A {2) E L - X2 - 1 
+b\E Xl A^ E X2 ~ Xl ~ x A^ E L ~ X2 ~ X . 



(27) 



(28) 



(29) 



(30) 



A, 



L-1 



E xi -1^4(1) e l ~ Xi ~ 1 A^ 



c 2 b 2 A^E Xl - 2 A^E L - Xl 
+ Cl c 2 A {1) E Xl - 2 A {2) E L - Xl 
+b 2 b 1 A < - 2 '>E Xl - 2 A^E L - Xl 
+c 1 b 1 A < - 2 '>E Xl - 2 A^E L ~ Xl , 



(31) 



A, 



L-2 



E xi-l A (2) E L-x x -l A {2) 



^1)^,-1 A {l) E L-x^-l 

+b 1 c 2 A^E Xl - 1 A^E L - Xl ~ 1 
+c 2 6i^ (2) E* 1 - 1 A^ E L - Xl ~ x 
+b\A^E^- 1 A^E L - x ^ 1 . 



(32) 



On the other hand, the configurations \x\, 2; x 2 = x\ + 
t + 1, 1 >, where we have collision of the arrows, give us, 
instead, the equations 



A2 .E x ^- 1 A^E t A^E L ~ xl ~ t ' 1 



HE X1 A^ E 1 - 1 A^ E 



L—xi —t—1 



(33) 



where we have used the interaction energy e/ = 
given by . As in the case n = 1 we identify the matri- 
ces A a (a = 1, 2) as composed by two spectral dependent 
matrices A kl and A k2 : 



A^=Y.^A kt E^ 2 \ 



(34) 



where <p l a {oi,i — 1,2) are unknown scalars and the spec- 
tral dependent matrices obey the commutation relation 
with the matrix E 



EA k = e j A ki E, j — 1.2. 



(35) 



The relations l|34|l and l|35|l . when inserted in (|14fl with 
n = 2, relate the momentum of the eigenstate with the 
spectral parameters: 



p = ki+k 2 . 



(36) 



Using (|34fl - l|36[l in l|27|) - (|32|l the non-zero norm condition 
for the eigenstate gives, for any x\ and x 2 , a single set of 
equations for the scalars <j> % a (a, i = 1, 2): 
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b\4>\ 

t>2<t>2 



b 2 C2e lk * 
C2b 2 e ikl 



&2C1 
b 2 bie ik * 
c 2 cie lfcl 



C1&2 
cic 2 e lfc2 
6i6 2 e ifel 



ci&ie ife2 
b lCl e' kl 



c 2 e i(k 1+ k 2) C2fcie i(fei+fc 2 ) j^cae'Cfci+fca) b 2 e ;( fcl+fe2 ) 

(37) 



'102 
'2 01 
■^1 



This last equation can be rewritten in the tensorial 
form as 



A 2 / 

L-2 \ 
X 






'<; 


^1 




<8> 


.< 




" b 2 




Cl 






cie ik 


1 b 




L 


(8 



c 2 e 



6 2 

ik 2 



Cl 



)( 



'I 

'2 



(38) 



We can always choose one of the fugacities to be unity. 
If we set an = 1 we see from <|25[) that the eigenvalue 
A 2 (ki, k 2 ) is given by 



A 2 (fci,fc 2 ) = Ai(fci)Ai(fc 2 ), 



(39) 



where Ai(fc) is the one-particle energy (|26fl. and the 
scalars 0^, satisfy 



4 fei A fei =s(kj,h)A ki A kj , A 2 ki =0, (* = 1,2) (44) 
where 

s(/cj , fcj) = 

Ai(fci)Ai(fej)6i - Ai(fcj)(fe 2 fci - c 2 ci + ai) + ai6 2 
Ai(Aj)Ai(fej)6i - Ai(fcj)(6 2 6i - c 2 ci + ai) + ai& 2 ' 

(45) 

The spectral parameters k\ and fc 2 , which were free up 
to now, are going to be fixed by the boundary relation 
1)150. Inserting (|34[l in (|15|1 and using the commutation 
relations l|35|) wc obtain 

3,1=1 

e~ v E ^ a ^L 1 e tt,(ei - at) A fcj At,^ £ - 4 '. (46) 
3,1=1 

Making the permutation j <-> I on the right-hand side of 
the last equation and using 1441) we obtain the spectral 
parameter equations that fix k\ and k 2 : 



Aki L 



= [-=-) s(k 3 ,h), 3 + 1 = 1,2. (47) 



cm 



Ai(fci) - 6 2 



i = 1,2. 



(40) 



Up to now the commutations relations of the matrices 
A ki among themselves are not known. They are going to 
be fixed by the "colliding" equations (|33[) . Inserting 134|) 
111 and using H35|) we obtain 



The solutions (ki,k 2 ) of these equations when inserted 
in give us the eigenvalues A 2 (fci, k 2 ). 

n > 2. The general case follows straightforwardly from 
the n = 2 case. The equations with no collisions lead 
us to identify the matrices A^ of the ansatz (|12|) as 
being composed by n spectral dependent matrices A kj 
(j = 1, . . . ,n), 



2 

E 

i,3=l,^3 



E Xl - 1 A kl E- t A kj E L - Xl - 3t+1 = 0, 



(41) 



where A 2 is given by (|39|) . This last equation can be 
rewritten conveniently by the use of the relations 



A?(fe/) + &2Ai(fcj) = -Ai{kj)bie ik * +e k i(b 2 b 1 - c 2 ci)(42) 



with j = 1, 2, and 



A a =Y J ^A k] E 1 - 2 \ a -1,2, 
satisfying the algebraic relation 



(48) 



EA kj =e lk >A kj E. (49) 

The equation l|14|) relates the momentum p of the eigen- 
state and the spectral parameters: 



b\<j>l b 2 -A 1 {k 1 ) 



Ax{k 2 



(43) 



that follows from (|25|l and l|4(J|) . respectively. These last 
two relations imply the commutation relations among the 
spectral matrices 



P = J2 k 3- 
3=1 



(50) 



The equations comming from the configurations where 
two arrows are at colliding positions fix the commutations 
relations among the matrices A k . : 
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A kj A kl = s{kj,ki)A kl A kj , 



Al. 



0, j = l, 



,n,(51) 



where s(kj,ki) is given by (|45|l . The equations with more 
than two arrows at collisions give relations that are auto- 
matically satisfied, due to the associativity of the algebra 
of the matrices A kj (J = 1, . . . , n). The eigenvalues of the 
transfer matrix are given by 



A„(A: 1 ,...,fc n ) = JjAi(%). 



(52) 



The spectral parameters {k{\ are fixed by the boundary 
relation i|15|) and are given by the solution of the nonlin- 
ear set of equations 



TT fe ikj 



i=i 



j = l,...,n, (53) 



with s(kj,ki) given by (05j). This last equation repro- 
duces for t = the spectral parameter equations for the 
six-vertex model obtained in [jj by using the coordinate 
Bethe ansatz. 

Since p = Sj=i kj we can rewrite the spectral pa- 
rameter equations l|53[) . in the sector with a number 
n (n — 1,2,...) of arrows and momentum p = -^1- 
(i = 0,1, ...,L- 1), as 



D ikj (L-2tn) -2ipt 



Y[s(kj,ki). 



(54) 



Consequently the eigenvalues belonging to the sector la- 
belled by (n,p) of Td-_d °f the interacting five- vertex 
model with a parameter t (t — 0, 1, 2, . . .) are related to 
those of the standard six- vertex model (t = 0). The six- 
vertex model is defined on a lattice size L = L — 2nt 
and with a momentum dependent seam |13| along the 
vertical direction of Fig.^3. The same phenomenon also 
happens on quantum Hamiltonians with hard-core exclu- 
sion effects [l2| . 

We should mention that is also possible to reobtain 
our solution of the interacting vertex models by using 
the coordinate Bethe ansatz. The eigenfunctions we de- 
rived through the matrix product ansatz give the edu- 
cated guess that we should use for the coordinate Bethe 
ansatz. The spectral parameter equations we obtain are 
the same as (JS3}. 

Before closing this section let us give a possible repre- 
sentation for the matrices A and E of the ansatz (|12fl . 
For a given solution {k%, . . . , k n } of the spectral param- 
eter equations <|54|) . in the sector with n particles, the 
matrices E and {A kl , ■ • ■ , A kri } have the following finite- 
dimensional representation 



E 





-ikl 



A ki = 



1=3 + 1 



j-1 

® 

Lz=i 

1 
1 



s(kj,ki) 
1 



1 




(55) 



where s(kj, ki) are given by (|45|l and A is obtained from 
(|48|l . The dimension of the representation is 2" and the 
products appearing on the ansatz are traceless. If we 
want a formulation of the matrix product ansatz where 
the trace operation is used, as in the formulation H13|) . 
it is quite simple to produce the matrix fl p that gives a 
non-zero value for the trace. We should stress that in 
the original formulation of the ansatz |(| , it was required 
unnecessary additional algebraic relations of the matrices 
E and A with the boundary matrix f2 p . The representa- 
tions, in this case, are probably infinite dimensional. 



V. ROOTS OF THE SPECTRAL PARAMETER 
EQUATIONS 

In order to complete the solution of any integrable 
model we need to find the roots of the associated spectral 
parameter equations (Eq. (|53|l in our case) . The solution 
of those equations is in general a quite difficult problem 
for finite L. However numerical analysis on small lattices 
allows us to conjecture, for each problem, the particular 
distributions of roots that correspond to the most im- 
portant eigenvalues in the bulk limit (L — * oo). Those 
are the eigenvalues with larger real part in the case of 
transfer matrix calculations. The equations we obtained 
in the last section were never analyzed previously either 
for finite or infinite values of L. Even in the simplest case 
t = 0, where the model reduces to the six-vertex model, 
the spectral parameter equations obtained in , through 
the Bethe ansatz, were not analyzed in the literature. 
A detailed numerical analysis of the Bethe ansatz equa- 
tions for the row-to-row transfer matrix of the six- vertex 
model, or the XXZ chain was previously done [Til Il8| . 
However, as we saw in last section, the corresponding 
equations for the diagonal-to-diagonal transfer matrix of 
the six- vertex model arc quite different and a numerical 
study of these equations is still missing. 

In our general solution of last section we have, for ar- 
bitrary values of the interacting range t (t = 0, 1, 2, . . .), 
five free parameters: (a±, b\, 62, c-x, 02)- The particular 
case where we have no interactions along the diagonals 
(ai = 0) is special and is not going to be considered here 
(see [l5j for a dicussion of the parameterization in this 
case). 

In order to simplify our analysis we are going hereafter 
to restrict ourselves to a symmetric version of our model 
with only three free parameters (5, b, c) , namely, 



Oo = l, ai = S , 61 



6(5, C1—C2— cS. (56) 



9 



The parameter 6 gives, in the case where t = 0, the 
contribution to the fugacity due to an electric field on the 
symmetric six- vertex model. For general values oit,p = 
— In 5, plays the role of a chemical potential controlling 
the number of arrows in the thermodynamic limit. 

Instead of writing the spectral parameter equations in 
terms of the spectral parameters (ki, . . . , k n ) as in (15411 . 
it is more convenient to write these equations in terms of 
the variables Xj = Al -p - , with Ai(fcj) given by l|26[) . In 
this case the eigenvalues of Td-d are given by 



A n — 5 n \\ ■ ■ ■ A n , 



(57) 



where {Xj} satisfy 



Xjib-Xj) 



b(b-X 3 )~c< 



L-2tr 



(-) n+i n 



t yj XiXj - 2AXj + 1 
A/A,- - 2AA; + 1 



j = l,...,n, (58) 



;=i 



and we have introduced the anisotropy parameter 



2b 



(59) 



We see, from (|57H - <|59|) . that we have now only two 
free parameters: b and c. The interacting parameter 5, 
that gives the contribution due to the interactions among 
the vertices along the diagonal, does not appear on the 
equations l|58|l and 159fl . It only gives an overall scale for 
the eigenvalues as shown in (|57|l . Inspired on the usual 
parameterization of the six-vertex model [2| we express 
the parameters b and c in terms of the parameters a and 
7: 



These equations are quite distinct from the correspond- 
ing spectral parameter equations derived for the row-to- 
row transfer matrix of the six- vertex model 0, 0] ■ Since 
no numerical analysis of the roots for this type of equa- 
tions were reported in the literature, we made an exten- 
sive numerical study of them for finite values of L and 
several values of the anisotropy A. In the particular case 
where A = (7 = -|) these equations can be solved ana- 
lytically. Solutions of these equations are obtained by the 
Newton method by using the distribution of roots {<Ji} at 
A = as the starting point to obtain the corresponding 
roots at other values of A ^ 0. Our numerical analysis 
shows that the eigenspectrum of Td-d is formed by real 
or complex-conjugated pairs of roots ensuring that the 
partition function is a real number. We verified that 
the eigenvalue with highest real part, belonging to the 
sector with n arrows, is real and corresponds to a zero 
momentum eigenstate {p — 0). The distribution of roots 
{er} corresponding to this eigenvalue has a fixed imagi- 
nary part, that depends on a and 7, and a symmetrically 
distributed real part, i. e, 



Imfo) 
J = ^ ■ 



7 — a 



, Real(crj) = Re&\((T n ^j), 



(63) 



We have also verified, for all sectors, the occurrence of 
several other real eigenvalues. In these cases the corre- 
sponding roots {<Ji} have imaginary parts given either by 

2 ul 2 27 

In the next section we are going to explore the topology 
of the above distribution of roots in order to extract the 
bulk limit (L — > 00) of the interacting five- vertex model. 



VI. THERMODYNAMIC LIMIT 



smer sin 7 

= 0(7, a) = -J-J- — , c = c(7,cr) = 



sin(7 — a) 



A = — cos 7. 



sin(7 — a) ' 
(60) 



The fact that A is a real number implies that 7 is real 
for — 1 < A < 1, and pure imaginary for |A| > 1. Since 
the right-hand side of (|58|l is the same for all values of the 
parameter t it is interesting, as in the six-vertex model 
Q , to make the change of variables Xj — > Oj , where, 



sinh(i7 — <Tj 
sinh a-j 



j = 1, . . . ,n. 



(61) 



In terms of these new variables {<Xj} the spectral param- 
eter equations become 



sinh(i7 — <jj) sinh(i7 — ia — <Jj) 
sinh(crj) sinh(i7 + Oj) 

-A- sinh(cr 3 -<ri+ ij) . _ 
sinh(<Tj - ui - ij) ' 



L — 2tn 



-Ivpt 



(62) 



We shall obtain in this section the free energy f(j, a, S) 
of the interacting five-vertex model in the bulk limit 
L — > 00. The spectral parameter equations (|54(l or 
()62J) tell us that this interacting model with parameter 
t (t = 0, 1, 2, . . .), on a lattice of size L and density of 
arrows p = ^, is equivalent to a six- vertex model (t = 0) 
on an effective lattice size L' = L — 2nt and density 
p' = ji = L ™2tn • This means that in principle we can 
derive the thermodynamic properties of the interacting 
five-vertex models by considering the many sectors, with 
fixed density of arrows, of the six-vertex (t = 0) model 
defined in the helical geometry of Fig. 

Since the relationship among the interacting five- 
vertex models with distinct values of the parameters t 
occurs on the eigensectors with distinct numbers n of 
arrows, we need to evaluate the highest eigenvalue of 
Td-d on the eigensectors with arbitrary density of ar- 
rows p = j-. The free energy per site in the bulk limit 
/(7, cr, 5) will be obtained from the largest eigenvalue of 
Td-d. Due to (|57|l the eigensector containing this eigen- 
value will have a density of arrows p = J that depends 
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on the particular value of the parameter 5, i. e, 



f(l, a, S) = min p / p (7, a, 5), 



where 



f P = 



lim 

n — >oc\L — >oo,p— j- 



\3=1 



and {Xj} are the roots of the equations (|58|l correspond- 
ing to the largest eigenvalue of Td-d in the sector with 
n = pL arrows. 

The numerical analysis of previous sections (see Ij61|) 
and iJSsSJl) indicates that the roots {Xj} in [|(t5|| correspond 
to a state with zero momentum and are given by 



2L — n(2t + 1) of empty places excluded by the hard- 
core interactions of the arrows (holes) satisfies rih > n 
and rih < n < 2 \+\ > or equivalently, < n < or 



(()4) JTI < n < 2^ . If 2iTT < « < 2m, we start again 
our matrix product ansatz l]12[| by associating now the 
matrices and not to the arrows, as before, but 
to the empty places excluded by the hard-core sizes of 
the arrows (holes) . The spectral parameter equations we 
obtain are the same as ijo^)) and (|59")l . except that the 
number n of variables {<7j} and {vj}, on the right-hand 
side of these equations, is replaced by rih = 2L — n(2t+l). 

In order to calculate the free energy l(55)l. as usual, we 
assume that as L — * oo the roots {vj} tend to an uniform 
distribution on a symmetrical interval of the real axis, 
-Q < v < Q. Let us denoted by 2(L - 2nt)R(Q\v)dv 
the number of roots {vj} on the interval dv, in the sector 
with density of roots f, i. e, 



a , = v j 



sinh(i7 — tjj) 
sinh aj 
.7 — <J 



Vn-jJ = 1,. ..,71, (66) 



where Vj are real numbers. Inserting the roots (|66|l in 
(|62|l we verify that both sides of this last equation are 
unimodular. Using the relations 



R{Q\v)dv 



< p < 



2(L-n(s-l)) 2(l-pt)' - P - (*+l) ' 

2L-(2s-l)rt _ 2-(2t + l)p 1 < 2 I ' L > 

2(i-n(s-l)) — 2(l-pt) ' (t+1) <~ P — (2t+l) ' 

In terms of this density of roots the function / p in l|65|) 
is given by 



sinh(z — iw) 
sinh(^ + iw) 



0(z, to) = arctan 



2 tanh z cot w 



1 — tanh z cot 2 w 
we can replace (|5^1l by the equations 

{L-2nt){e x {v j ) + e 2 {v j )) = 

n 

2TTl j +J29 3 (vi-v j ), 3 = 1 



(67) 



(68) 



)**•)' "J 



(69) 



z=i 



where 



?i(«) = %,^), e 2 (i;) = e(«,^- £ : 



= 6(v,j), 



(70) 



and /j are integers or half-odd integers depending if n 
is odd or even, respectively. Our numerical analysis on 
small lattices indicates that those integers are given by 



h = i 



,n. 



(71) 



In deriving the equations l|62|) and l|69|l it is better to 
consider separately the cases where the number rih = 



f p = -2k B TplnS- 2k B T{\ - 2pt) 



Q 

x / R(Q\v)1n 
l-Q 



sinh(i 



(7+"-) 



V) 



sinh(i 



(7-gj 



+ v) 



dv. 



(73) 



In order to obtain R(Q\v) we notice from l|69|) and (|71|) 
that for a given value Vj of the variable Vj the integer 
Ij + gives the number of variables vi with I < j. 
Using this fact l|69|l can be written as 



(L - 2tn)(9 1 (v) + 6 2 (v)) = -7r(n + 1) 

+4tt(L - 2tn) I R(Q\v')dv' 
J-Q 

+2{L-2tn) / 9 3 {v - v)R{Q\v')dv '. (74) 
J-Q 

Taking the derivative of this last expression, with respect 
to v, we obtain equations that are conveniently expressed 
in the trigonometric regime, —1 < A = — cos 7 < 1, by 



R(Q\v) = 
1 sinfer 



7) 



sin(7 



2-7T cosh(2w) — cos(er + 7) 
1 f Q sin(2 7 ) 



cosh(2(w — v')) — cos(27) 



2-7T cosh(2w) — cos(7 — a) 
R(Q\v')dv'. (75) 
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In the hyperbolic regime, where A < —1, it is better to 
change 7 — > — iX, a — > — icr, 1; — * —iv, and the equations 
are 



1 



sinh(A — a) 



R(Q\v) = 
1 sinh(cr + A) 

2n cosh(cr + A) — cos(2w) 2tt cosh(A — a) — cos(2i;) 
1 f Q sinh(2A) 



7r J_q cosh(2A) — cos(2(v — v')) 
A = - cosh A < -1. 



jrrR(Q\v)dv', 



(76) 



The coupled integral equations lfT2l and (|?5jl or ff^ 

and l|76[l give us, for a given density of roots p, the param- 
eter Q = Q(p) and the density of roots R(Q\v). Solving 
numerically l|72|) and (|75|l or (|72|l and l|76|) and inserting 
the result in i|73fl . we obtain f p — f p {p/,a,5). It is im- 
portant to stress (see i|64|l) that for given values of the 
parameters 7, a and 5 the free energy is given by the 
minimum value of f p . In Fig. [7| we show, for 7 = § 
(A = — |), cr = 0.1, 6 = 1 the curves f p for some val- 
ues of the parameters t. As we see on this figure the 
minimum value of f p that gives the free energy occurs, 
for fixed values of the parameters 7, a and d, on sectors 
whose density of arrows p depends on the parameters t 
(t = 0, 1, 2, . . .). In the case where t — the choice 5 = 1 
gives us the symmetric six- vertex model on the helical ge- 
ometry (Fig.^>). The largest eigenvalue of Td-d occurs, 
in this case, in the sector with density p = 1. In Fig. [S] 
we show, for some values oft, the free energy / = min p / p 
as a function of the parameter 5. 




(2t+1)p 



FIG. 7: Functions fp/KsT, defined on l|65^ . as a function 
of p, for the interacting vertex model with t = 0, 2, and 4. 
These curves are evaluated for A = — |, a = 0.1 and 5 = 1. 
In this figure the values on the horizontal axis are multiplied 
by (2t + 1) to represent the curves on the same scale. The 
minimum values of these functions give the free energy. 

As we have already mentioned the spectral parameter 
equations (|58|l were not analyzed even for the standard 
symmetric six-vertex model. This case corresponds to 
the case where 6=1 and t = (s = 1). Our numerical 
results (see e.g. Fig. [7J indicate that in this case the 
parameter Q in (JJ^J-lJJSJl takes the value Q — > 00 and 
Q = 7T for the trigonometric (—1 < A < 1) and the 



K„T 



1 

t=4 


V c ' : -v.:/. 1~ 


- \ t=2 




\ t=0 







1 



FIG. 8: Free energy, in units of KbT, as a function of 8, 
for 5 > 1 for the interacting vertex models with t = 0,2, and 
4. The curves shown are for the anisotropy A = — i and 
parameter a = 0.1. 



hyperbolic (A < —1) regimes, respectively. We can use in 
this case Fourier transforms in order to solve the coupled 
integral equations and GSJ or lO and (J7HJI. After 
some straightforward algebra we obtain for — 1 < A = 
— COS7 < 1, —7 < a < 7, 



R(Q\v) = R(oo\v) 



n cosh(^)cos(^) 
4 7 cosh(^) + cos(^)' ( j 



that gives the free energy 



/ = ~k B T 



2 cosh(^)cos(^) 
7cosh(^) + cos(^) 



x In 



cosh(2'y) — cos(7 + a) 



dv. 



(78) 



On the other hand, for A = 
-A < a < A, 



cosh(2w) — cos(7 — cr) _ 

coshA < -1, A > 0, 



R(Q\v 
1 



R(ir\v) 



coshf — Isech — )e 



'dx 



and 



f = -k B T(a+J2 



(-iy 



sinh(2TOcr)sech(mA)e mX ). 



In table we present on the last line the numeri- 
cal evaluation of the integral H78|) for some values of 
A = — cos 7 < 1 . In order to compare with the re- 
sults for finite values of L we also show on this table the 
results obtained by solving directly the spectral parame- 
ter equations 1)62(1. The free energy per site in the bulk 
limit should be independ of the geometry of the lattice. 
Indeed the results of table [I] obtained for the symmetric 
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six-vertex model in the helical geometry, coincide with 
those obtained from the row-to-row transfer matrix in 
the usual geometry (Fig.GJi). It is however interesting to 
observe that although the final results are the same, the 
integrand in (|78() for the helical case (Fig. QJj) is distinct 
from the corresponding one in the geometry of Fig. 
(see Eq. (8.8.19) in @) 



L 


7=¥ 


7=f 


7=1 


6 


0.03512818236 


0.06453944019 


0.11709717612 


10 


0.03477636310 


0.06406728711 


0.11638357858 


18 


0.03464007645 


0.06388499114 


0.11610862696 


34 


0.03459635046 


0.06382657413 


0.11602059391 


66 


0.03458384714 


0.06380987622 


0.11599543821 


130 


0.03458049496 


0.06380539993 


0.11598869526 


oo 


0.03457933094 


0.06380384561 


0.11598635395 



TABLE I: Free energy of the six-vertex model (t = 0) on the 
distorted lattice of Fig. for some values of the anisotropy 
A = — cos 7 and lattice sizes L. The values are obtained by 
solving 1621 with a — 0.1 and using 1611 and 16511 with 5 = 1. 
The last line gives the bulk limit obtained by solving (1781 . 



VII. CONFORMAL INVARIANCE AND 
OPERATOR CONTENT 



In this section, by exploring the conformal invariance 
of the infinite system, we shall obtain the critical behav- 
ior of the massless phases of the interacting five-vertex 
models introduced on this paper. The results of the last 
section tell us that for A < 1 the eigensector of Td-d 
presenting the largest eigenvalue, changes continuously 
with the value of the interacting parameter 6. This imply 
that for A < 1 the introduced interacting models with 
arbitrary values ofi = 0,l,2,... are on a critical phase, 
as we change continuosly the parameter S in the bulk 
limit (L — > oo). This behavior for the case t = (s = 1) 
is expected since, as we discussed in earlier sections, the 
interacting model recovers the six-vertex model. In this 
case the interacting parameter 6 plays the role of a chem- 
ical potential (or electric field) controlling the density of 
arrows. The critical behavior of this six- vertex model fol- 
lows from that of the XXZ chain in the presence of an 
external magnetic field 16]. The region where A > 1 
(ferroelectric region) is non-critical for any value of t and 
5. The largest eigenvalue of Td-d belongs to the sec- 
tor with no arrows (p = 0), for S < 1, or to the sector 
with the largest possible density of arrows (p = ^rpr ) f° r 
S> 1. 

The conformal anomaly and anomalous dimensions 
(related to the critical exponents) of the underlying con- 
formal field theory governing the critical fluctuations on 
the critical phase of our models are going to be calcu- 
lated from the finite-size corrections of the eigenvalues 
of the transfer matrix Tu—Di or the related Hamiltonian 



H = — ]xi(Ti}-d)- The leading behavior, as L — > oo, of 
the smallest eigenvalue F (L) of H (ground-state energy) 
El: 



F (L) =F {oo) 



irv s c 
~~6T 



o{L~ 



(79) 



gives the conformal anomaly c and the sound velocity v s . 
Moreover, to each primary operator <f> with dimensions 
x<fr, in the operator algebra of the critical system, there 
exists an infinite tower of eigenstates of H = — ln(2/j_£>), 
whose eigenvalues F^ , (L) behaves asymptotically as 

El 



fLw(L) = F (L) + ^(^ + m + m') + o(L- 1 ),(80) 

with m, m! = 0,1,2.... 

Before considering our model with general values of 
t (t = 0, 1, 2 . . .) and 6, let us consider the special case 
where t = and 5 = 1. In this case we recover the 
standard symmetric six-vertex model in the absence of 
electric fields, but for the helical geometry of Fig. ^d. 
Similarly, as in the standard geometry of Fig.^, |21[. our 
numerical and analytical analysis indicate that the model 
is critical for — 1 < A < 1, and massive for |A| > 1. In 
agreement with earlier studies of the six-vertex on the 
standard geometry, and of the XXZ quantum chain |14| , 
the critical fluctuations in the region — 1 < A = — cos 7 < 
1, are governed by a conformal field theory with central 
charge c = 1 and anomalous dimensions given by 



= n 2 x p + - — , n, m = 0, ±1,±2, (81) 

4Xp 



where 



7T — 7 

2tt 



(82) 



Our numerical and analytical analyses show that dis- 
tinctly from the row-to-row transfer matrices in the 
standard geometry where the sound velocity is unity 
IS l2fl |. in the helical geometry one has v s = 
u s (7, a) ^ 1. The sound velocity from (fTO - )) can be calcu- 
lated from the bulk limit of the finite-size sequence 



[Fo(D 



F(oc)]^. 

7T 



(83) 



For A = (7 = 77) we can compute analytically these 



finite-size corrections and obtain 



1 a 
— tan — . 

2 2 



(84) 



For A / our calculations were done numerically and 
we do not have an analytical form for v s . In table ITT1 we 
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show some of the estimates v s (L) given by (|83[) obtained 
for some values of A at a = 0.1. In order to illustrate 
our numerical evaluation of the dimensions 1)81(1 we show 
in table Hill for a = 0.1 and some values of 7, the finite- 
size estimates for x p , together with the expected results 
(JH2J). These estimates are obtained by using in lj%0"jl the 
largest eigenvalues of Td-d in the sectors with n = L 
and n = L — 1 arrows. 



where, for — 1 < A = — cos 7 < 1, 



L 


7=¥ 


7=f 


7=f 


6 


0.03773624455 


0.05057575759 


0.07637450614 


10 


0.03763037006 


0.05031361905 


0.07586431519 


18 


0.03758898261 


0.05021239628 


0.07566185294 


34 


0.03757566153 


0.05017995742 


0.07559488486 


66 


0.03757184873 


0.05017068491 


0.07557509072 


130 


0.03757082640 


0.05016819918 


0.07556960190 



TABLE II: Finite-size estimators v s (L) for the sound velocity 
of the six-vertex model (t = 0) on the distorted lattice of 
Fig-Q 3 - The values are given for some values of the anisotropy 
A = — cos 7 and lattice sizes L. 



L 


7=¥ 


7=f 


7=f 


6 


0.16676687232 


0.24712196038 


0.32156619180 


10 


0.16669813477 


0.24895867881 


0.32848712047 


18 


0.16667118823 


0.24967533316 


0.33163455071 


34 


0.16666377768 


0.24990588089 


0.33279409868 


66 


0.16666111003 


0.24997186077 


0.33308163424 


130 


0.16666077565 


0.24998955436 


0.33328029440 


Xp 


0.1666.. 


0.25 


0.3333.. 



TABLE III: Finite-size estimators for the critical exponent 
x p of the six-vertex model on the distorted lattice (t = 0) of 
Fig- 03- The values are given for some values of the anisotropy 
A = — cos 7 and lattice sizes L. These estimators are obtained 
by solving 162H . 

Let us consider the general case where t — 0, 1, 2 . . . 
and 5 > 0. In this case we have a massless phase for 
A < 1. For a given value of S the density of arrows p, 
in the bulk limit, is obtained by the condition = 0, 
where f p is given by (|73|) . As we change <5, from zero to 
infinite, p increases from to 2 t+i • ® nr numerical anal- 
ysis indicates that the underlying conformal field theory, 
as in the case t = and S = 1, has also a central charge 
c = 1 and the Coulomb gas type of conformal dimensions 
given in l|81(l , but with a value of x p that depends on (5,7 
and a. The parameter x p , that can be estimated from 
the finite-size corrections of the smallest eigenvalues of 
H = — ln(Trj)— d)> can be calculated analytically. This is 
done by applying to our relations lf?2" )l -l(?rj |l of last sec- 
tion, the method used in for the XXZ in a magnetic 
field. We obtain 



V (Q\v) = l- 



1 



Q 



sin(2j)r)(Q\v') 



7r J_q cosh(2(u — v') — cos(27) 
and, for A = — cosh A < —1, 



dv', (86) 



V (Q\v) = l-- 

7T 



9 sinh(2A)?7(Q|?/) 
q cosh(2A) — cos(2(t> — v') 



-dv'. (87) 



For a given density of arrows p the parameter Q — Q(p) is 
obtained by solving the coupled integral equations l|72|) 
and J7SJ) for -1 < A < 1, and and J7U for A = 
— cosh A < — 1. Similar integral equations ha ppe n for 
the XXZ chain with hard-core exclusion effects [231 ]. 

At A = Eq. dHSJ gives rj(Q\v) = 1, and from (gZJ we 
obtain the curve 



1 



4(1 -tpf 



0. 



For fixed values of t, a and 7 (or A) the interacting 
parameter S will fix the density of arrows p as well as the 
parameters Q = Q(p) appearing in (|72|) . l (75)l (or l(7B)ll 
and l|85() - (|87|l . In Fig. [!|]we show the curves of Xp as a 
function of p at A = — h and for some fixed values of 
t. In Fig. ITTil and Fig. ^| the values of function 
of p are shown for the interacting model with parameter 
t = 2. The curves of Fig. ED an d Fig. are for the 
regimes — 1 < A < 1 and A < — 1, respectively. 

As we see in Fig.llOland Fig.^Jthe curves have a com- 
mon intercept x p = 1/4 at p — and x p = 25/4 = 6.25 
at p = 2/5 = 0.4. In fact our numerical results, for arbi- 
trary values oft and A, shows that at low density (p — ► 0) 
and high density of arrows (p — * 2 t+\ ) ^ ne anomalous di- 
mension x p take the values x p = 1/4 and x p = ^"t 1 -* , 
respectively. This happens since at those extreme lim- 
its either the density of arrows or the density of holes is 
low and the anisotropy A, that is only important for the 
arrows (holes) at the closest positions, plays no role. 

We also see in Fig. ^] the appearance of a peak 
(x p = 6 2 /4 = 9) at the density p — \, that is independent 
of the anisotropy A. For general values of t and A < — 1 
this peak occurs at the " half- filled" density p = for 
any A < — 1. This happens because, for A < —1, the 
system at this " half- filling" has a mass gap. An infinites- 
imal change of p (induced by a change of 6) destroys this 
mass gap. Since around this density the main effect is 
the disappearance of the gap, the critical exponent x p 
does not depend on the A value. A similar effect is also 
expected in the region A < — 1 of the XXZ or six- vertex 
model on the presence of an external magnetic field. 



VIII. CONCLUSIONS 



1 



(1 



t P )- 2 r 2 {Q\v) 



We have introduced in this paper a special family of 

(85) 

exact solvable interacting five- vertex models. Besides the 
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2.5 - 




(2t+l)p 

FIG. 9: The critical exponent x p , as a function of p, for 
the interacting vertex models with A = — | and interacting 
parameter t — 0,2 and 4. In the figure the horizontal and 
(vertical) axis are multiplied (divided) by 2t + 1 in order to 
represent the curves on the same scale. 




FIG. 10: The critical exponent x p , as a function of p, for 
the vertex model with interacting parameter t — 2 and some 
values of the anisotropy A = — cos 7 < — 1. 




usual nearest-neighbor interactions imposed by the lat- 
tice connectivity, such models also present hard-core in- 
teractions along one of the diagonals of the lattice. The 
range of the additional interactions depend on a fixed pa- 
rameter t (t = 1, 2, . . .). These vertex models can also be 
interpreted as if the arrows entering in the vertex con- 
figurations have an effective hard-core size s = 2t + 1. 
This interpretation allows us to identify the extension of 
these models to the case t = (s = 1) as the standard 
six-vertex model. In this sense, although our interacting 
models have only five vertex, they are generalizations of 
the standard six-vertex model. 

Although the exact solution of the transfer matrix are 
usually obtained through the Bethe ansatz, we have de- 
rived the solution of these models by using the matrix 
product ansatz introduced in [||, to perform the exact 
diagonalization of quantum chains. This first applica- 
tion indicates that this ansatz is also suitable for transfer 
matrix calculations. 

Due to the extra interactions along the diagonals of the 
lattice, the solution of these models was simplified by con- 
sidering the diagonal-to-diagonal transfer matrix. A nu- 
merical analysis of the spectral parameter equations for 
small lattices enabled the calculation of the free energy of 
the models in the bulk limit. This calculation was done 
on a symmetric version of the model that recovers for 
t = the symmetric six- vertex model. Our results shows 
that the models, for any value of t (t = 0, 1,2,...), exhibit 
a massless phase for arbitrary values of the anisotropy A 
given by l|59|l . Exploring the consequences of the con- 
formal invariance of the infinite system we also calcu- 
late, on the massless phase, the conformal anomaly, the 
sound velocity and the operator content of the underlying 
field theory governing the long-distance critical fluctua- 
tions. Underlying these models there is a c = 1 conformal 
field theory whose operators have a Coulomb gas type of 
anomalous dimensions x n ^ m (see (|81fl ). The compactifi- 
cation ratio defining these dimensions varies continuosly 
with the free parameters of the model. Our analytical 
and numerical calculations show that the effect of the in- 
teracting parameter t is to increase the critical exponents 
as t increases, while the extra interaction energy, related 
to the parameter 5, plays the role of a chemical poten- 
tial controlling the number of arrows on the ground state 
eigenfunction. 
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